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Abstract 

Incompressible stationary flows of ideal plasma are observed. By introduc- 
tion of curvilinear system of coordinates in which streamlines and magnetic 
force lines form a family of coordinate surfaces, MHD equations are partially 
integrated and brought to a certain convenient form. It is demonstrated 
that the admissible group of Bogoyavlenskij's symmetry transformations [1] 
performs as a scaling transformation for the curvilinear coordinates. An- 
al5d;ic description of stationary flows with constant total pressure is given. 
It is shown, that contact magnetic surfaces of such flows are translational 
surfaces, i.e. are swept out by translating one curve rigidly along another 
curve. Explicit examples of solutions with constant total pressure possessing 
a significant functional arbitrariness are given. 

Key words: ideal magnetohydrodynamics, stationary flows, exact 
solutions, contact magnetic surface, translational surface. 



Introduction 

Ideal magnetohydrodynamics equations (MHD) describe macroscopic mo- 
tions of infinitely conducting plasma under the action of internal pressure, 
magnetic and inertial forces. MHD approximation is applicable to wide 
classes of physical phenomenon all the way from plasma confinement to astro- 
physical problems [2, 3]. Well-studied in literature linear and low-dimensional 
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reductions usually have limited areas of application [4]. At the same time, 
the numerical analysis of MHD equations is complicated by the essential 
multi-dimensionality of MHD processes (e.g. the turbulent dynamo problem 
[5]), and presence of various types of strong and weak discontinuities [6, 7]. 
In this connection, analytical investigations based on construction and in- 
terpretation of exact solution to ideal MHD equations play significant role. 
Of course, it is impossible to construct (or even to define the notion of) the 
general solution to MHD equations. However, particular exact solutions give 
opportunity to describe main features and singularities of plasma flows on 
the level of explicit analytical relations. 

In this paper we observe the stationary ideal MHD equilibria equations: 

(v- V)v-(B- V)B + VP = 0, 

(B-V)v = (vV)B, (1) 

divv = 0, divB = 0. 

Here v is the velocity vector, B is the magnetic field, P = p + ||Bp is the 
total pressure, p is the hydrodynamical pressure of plasma. The dot denotes 
the scalar product in Euclidean space. Plasma density is assumed to be 
constant and is taken to be unity: p = 1. 

A significant part of known solutions to equations (1) belong to classes 
of either pure ideal fiuid (B = 0) [8] , or static equilibrium ( v = 0) [2] , or to 
the intermediate case of field-aligned flows (v = A(x)B) [9, 10, 11], which is 
known to be reducible to one of two previously mentioned [12]. Assuming 
the set of field-aligned fiows to be well-analyzed, we observe the general case 
of plasma fiows with non-collinear fields v and B. 

The induction equation (the second equation of (1)) states, that vector 
fields V and B commute. In sections 1, 2 we introduce a curvihnear sys- 
tem of coordinates, which reduces MHD equations (1) to a certain simpler 
form. A great deal of information about the fiow is provided by the con- 
tact magnetic surfaces weaved out of streamlines and magnetic lines of the 
fiow. Each of contact magnetic surfaces can be taken as an infinitely con- 
ducting impermeable wall bounding the fiow. In section 3 we demonstrate, 
that contact magnetic surfaces are constructed in the curvilinear coordinates 
automatically part of the solution of eqiiations. 

The remarkable property of equations (1) is that they have an infinite- 
dimensional symmetry group of Bogoyavlenskij's transformations [1, 12]: 

V = Bsinh/((^) -I- vcosh/((^), B = Bcosh/((^) -|- vsinh/((^). (2) 



2 



Here Lp is an arbitrary function, which conserves along streamhnes and mag- 
netic force hnes of the flow: 

(vV)v? = 0, (B-V)^ = 0. 

The group allows reproduction of new exact solutions from known ones [13] , 
construction of flows with current sheets and non-symmetrical equilibria 
states [12]. It is known [14] that Bogoyavlenskij's transformations preserve 
the stability properties of the solution. In section 2 we demonstrate that 
Bogoyavlenskij's transformations generate certain scale transformations for 
our curvilinear coordinates. Choice of the discontinuous scale factor allows 
producing a solution with a current sheet from any continuous solution of 
equations (1). 

Using the reduced system of MHD equations in section 4 we describe flows 
of incompressible plasma distinguished by the constant total pressure condi- 
tion: p+ ||Bp = const. Construction of these solutions requires integration 
of an overdetermined system of partial differential equations, obtained by 
neglecting VP term in equations (1). Note, that isobar ic (p = const) and 
barochronous (p = p{t)) solutions to pure gas dynamics equations (B = 0) 
are completely described in [15, 16, 17]. The method of integration used in 
these articles was based on the analysis of algebraic invariants for the Ja- 
coby matrix 9v/9x. In MHD equations there are two vector fields v and 
B, and, correspondingly, two Jacoby matrices. The compatibility conditions 
of the overdetermined system are formulated in terms of algebraic invariants 
of both matrices and of their products in various powers. This makes the 
analysis cumbersome and does not allow obtaining the closed form of the 
solution. 

In contrast to the described approach, in the curvilinear coordinates the 
class of solutions with constant total pressure is intrinsic and is described in 
section 4 by explicit formulae with significant functional arbitrariness. We 
prove, that every contact magnetic surface of stationary incompressible MHD 
fiow with constant total pressure is a translational surface, i.e., is obtained 
by a parallel shift of one curve along another curve in 3D space. Both curves 
are subjected to some restrictions. In section 5 we demonstrate examples 
of exact solutions describing flows in curvilinear channels. The functional 
arbitrariness of the solutions allows varying of the flow picture in wide extend. 
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1. Curvilinear coordinates 

In order to simplify the form of MHD equations let us introduce the 
following vector fields: 



or equivalently 
Equations (1) give 



B, b = v + B, (3) 



v=i(b + a), B = i(b-a). (4) 



(a- V)b+(b- V)a + 2VP = 0, (5) 
(a- V)b= (b- V)a, (6) 
diva = 0, divb = 0. (7) 

Equation (5) by virtue of (6) transforms to 

(a-V)b + VP = 0, (8) 

With new unknowns, transformation (2) takes the following form: 

a = (/?(x) a, b = h. 

Here (p is the arbitrary function, which conserves along vector fields a and b: 

(a-V)(/? = 0, (h-V)ip^O. 

Equation (6) is equivalent to the vanishing of the commutator of vector 
fields a and b: 

[a, b] := (a • V)b - (b • V)a = 0. (9) 

Under assumption of linear independence of these vector fields the commu- 
tativity condition (9) implies, that the fields can be taken as a coordinate 
basis, i.e. integral curves of the fields can serve as coordinate curves of a 
curvilinear coordinate system [18]. Indeed, let us observe curvilinear co- 
ordinates {k^,k'^,k^) such that k^- and /c^-coordinate curves coincide with 
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integral curves of vector fields a and b correspondingly. In other words, in 
the space let us define a difeomorphism 



x = x(k), det(^^^^0 (10) 

such that 

d'x. 

Compatibility of equations (11) for functions x(k) follows from equations (9). 
Indeed, calculation of the mixed derivative gives 

/9x\ 5 / / x\ 9a 9x 9a, 



a^x 


d 






92x 


d 


dk^dk"^ 


dk^ 



/9x\ 9 , . dh dx. dh 

[WJ^W (^^"^^ = d^dk^ = ^ = ■ 



(12) 



Right-hand sides of these equations coincide by virtue of equations (6), which 
imphes the compatibihty of equations (11). In variables x(k) the induction 
equation (6) is identically satisfied. 



2. Equations in the curvilinear coordinates 

In order to rewrite equation (8) in the curvilinear system of coordinates 
it is required to calculate vector VP in variables k. Note, that vector VP in 
fact should be treated as a covcctor, transformed to vector with the use of 
the identity metric tensor. In other words, column-vector VP in equations 
(8) must be correctly written as a transposed row-vector (covector): (VP)*. 
Let us denote by the lower index at symbol V the set of variables, to which 
the gradient is taken. We have 

By multiplication of equation (8) from the left on the transposed Jacoby 
matrix and taking into account formulae (12) we obtain 
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In the coordinate representation equation (13) have the form 

^WdkW^W~' ^-l,2,3. 

It is remain to transform equations (7) to curvilinear coordinates (10). To this 
end we use the invariance of differential 2-forms under any non-degenerate 
change of variables [18]. Let us observe the following 2-forms: 

cu^ = a^dx'^ A dx^ + a^dx^ A dx^ + a^dx^ A dx'^, 

a;2 = b^dx"^ A dx^ + b'^dx^ A dx^ + b^dx^ A dx'^. 

Here a* and 6* are coordinates of vectors a and b in the Cartesian frame of 
reference. Equations (7) are equivalent to 

rfo;^ = 0, du^ = 0. (14) 

By rewriting 2-forms uj^ and a;^ in variables (10), (11) we obtain 



Q dk^ Adk^, ^ Q dk^ Adk^, det ( ^ ) . 



V(9k J 

Closedness conditions (14) are equivalent to the following equations: 

Thus, the original system of equations (1) is reduced to equations (13), (15) 
for unknowns x = x(k), P — P(k), and to representations (4), (11) for 
sought vector fields v and B. Equations (15) are integrated as q = f{k^) 
with arbitrary function /. Note, that by virtue of invariance of equations 
(13), (15) under the change of variable — >• g{k^), by the suitable choice of 
function g one can always make f{k^) — 1. In the final form the resulting 
system of equations reads 

dx dP 

0, (16a) 
0, (16b) 
0, (16c) 



dk^ 


dk^dk'^ 




dk^ 


dx. 


d'^x 


+ 


dP 




dk^dk"^ 


dk"^ 


dx 


d'^x 


+ 


dP 



dk^ 



det ( ^ ) = 1. (16d) 



System (16) contains less number of equations and less number of unknowns 
than the original MHD equations (1). 

Note, that system of equations (16) is invariant under the change of pa- 
rameters k of the following form: 

~k'^^{k')k\ ~k' = -^, k'^k' (17) 

with arbitrary function (/?. This transformation follows from Bogoyavlcnskij's 
transformation (2), admitted by the original system of equations (1). From 
the point of view of introduced parametrization, formulae (17) specify a 
scaling transformation for parameters k^ and k'^. Besides, system (16) is 
invariant under transformations 

k^^k^+^{k^), k^ = k^ + x{k^) (18) 

with arbitrary functions i/j and x- In what follows we threat transformations 
(17), (18) as equivalence transformations on the set of solutions of equations 
(16). 



3. Geometrical interpretation 

Each dependence x = x(k), P = P(k) satisfying equations (16) deter- 
mines some incompressible flow of infinitely conducting plasma. A streamline 
of the flow X = Xsi and a magnetic force hue x = x^^, passing through point 
x(ko) are parametrically defined as 

^si = X ^/cq -|- -s, /cq + -s, k^ , = X ^/cq — —s, kg + —s, k^ , 

s e [0, -S] C R (19) 

Coordinate surfaces k^ — const are weaved out of streamlines and magnetic 
hues of the flow, hence, they can be treated as contact magnetic surfaces. 
Each of these surfaces can be taken as rigid infinitely conducting boundary, 
limiting the domain of the flow. 

Note, that coordinate surfaces k^ = const provide an explicit form of 
the fibration of Euchdian space by a family of surfaces ^(x) = const 
introduced in section 2 of paper [1]. This also allows us to give an exphcit 
description of flows with contact sheets presented in [1, 12]. Indeed, Bogoy- 
avlcnskij's transformation (17) does not change contact magnetic surfaces. 
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although it transforms vectors a and b according to relations (11), (17). 
The total pressure P is the invariant of the transformation. Thus, by acting 
with the transformation (17) on any continuous solution of equations (16) 
and taking function (p{k^) with jump discontinuity at some value = c, 
we obtain a new solution, which contains a contact discontinuity at the sur- 
face X = x(A;^, A;^, c). The surface of contact discontinuity carries the electric 
current with the density [19] 



Here n the unit normal vector to the discontinuity surface and directed to- 
wards side "2"; B^'^ are limiting values of magnetic field vector B on the 
discontinuity surface approached from side i. With the use of formulae (11) 
and transformation (17) we obtain the following formula for the surface cur- 
rent: 



Hereafter, the lower index i denotes the partial derivative with respect to k^; 
sign X stands for the vector product in M^. All vectors in the equation (20) 
are calculated at k^ = c; parameters k^ and /c^ are free. By (^^'^ we denote 
the limiting values of function (p{k^) at A;^ ^ c ± 0. 

4. Solution with constant total pressure 

Let us suppose, that total pressure P is constant in all space, occupied 
by the flow: 



J = nx(B(2)-B«). 




(20) 



P = const. 



In this case equations (16a)-(16c) are equivalent to 



dk^dk'^ 

which implies the following representation of vector x: 



(21) 



Substitution of the representation (21) into equation (16d) gives 

(Ti • (t2 X (0-3-^x3)) = 1. 



(22) 
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Thus, the description of the sohitions with constant total pressure is reduced 
to the separation of variables k^, in equation (22). 

Contact magnetic surfaces for solutions with constant total pressure are 
defined by the equation (21) at fixed value of k^. By virtue of this repre- 
sentation, these surfaces belong to a class of translational surfaces [20]. By 
definition, translational surfaces are those, which can be obtained by the 
parallel shift of one fixed curve in 3D space along another fixed curve. In the 
general form translational surfaces were analyzed in classical works by Lie, 
Darboux, Wirtinger, and Poincare (see [21] and citations therein). As it was 
first noticed by S. Lie [22], any translation surface, which could be swept out 
in more than one way by translating one curve rigidly along another curve 
is completely determined by a selection of a fourth-order algebraic curve on 
the plane. The cited works give classification of translational surfaces, which 
have several independent parametrizations of the form (21) with fixed value 
of k^. Translational surfaces also arise in study of billiards in rational poly- 
gons [? ]. To any such polygon there corresponds a unique translational 
surface such that the billiard fiow in the polygon is equivalent to a geodesic 
fiow on the surface. As far as we know, the relation between translation 
surfaces and fiows of infinitely conducting plasma was not yet discussed in 
literature. 

5. Examples of solutions 

The general case of separation of variables in equation (22) is cumber- 
some, therefore in this paper we limit ourselves to the case of special depen- 
dence (21), where vector field a does not depend on k^: 

x = (T(fc^) + t(A;^A;2). (23) 

In what follows, we denote by ei, e2, and 63 a triple of orthonormal constant 
vectors in M^. The decomposition of vector r in the basis of vectors have 
the following form: r = ^i=iT*ej. 

The dimension of the linear space {cr}, spanned by vectors cr(k^) for 
various values of /c^, can be 1, 2, or 3. Below we observe each of three cases 
separately. 

a) dim{cr} — 1. In this case vector cr can be represented as 

cr — a{k^)ei. 
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By plugging this representation into (22) and by using the decomposition of 
vector T we obtain 



J.2 3 

^2 ^2 

^3 ^3 



1. 



Hence, accurate to equivalence transformations (17), (18) this equation im- 
phes a{k^) — k}, and T2TI — t^t^ — 1. Finally, we obtain the following 
solution: 



x = fc^ei + T(A;^,A;^), 



T"2 ^ 
^1 ^1 



1. 



(24) 



The function T^[k'^,k^) remains arbitrary. 

b) dim{cr} = 2. We choose the following representation of vector a: 

(T^a{k')ei + p{k')e2 

with linearly independent functions a and (3. By plugging this representation 
into equation (22) we find the following relation: 



a'{k^) 



J.2 ^3 

^2 ^2 

T-2 3 

^3 ^3 



+ (5\k') 



T-3 1 

^2 ^2 

T-3 1 

^3 ^3 



= 1. 



By virtue of linear independence of functions a and P this relation is satisfied 
(accurate to the equivalence) only when a — k^, rfrg — T2t| = 0, and 
— 1- The second equation implies = I^i'T^) with arbitrary 



T-2 3 

'2 '3 



function F. Hence, we obtain the solution 

x= (^k' + F{T%k',k'))y, + (^P{k') +T\k^k')y2 + r'{k',k')es, 



T-2 3 

^3 ^3 



(25) 



Here (3 and F are arbitrary functions of their arguments, 
c) dim{<T} = 3. We use the following representation of cr: 

a = a{k^)ei + i3{k^)e2 + ^{k^)es 

with linearly independent functions a, P, and 7. Prom equation (22) we 
obtain 



a'{k') 





^2 




'2 












^3 


^3 




rl 





1. 
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Accurate to the equivalence, this implies a — k^, 




0. 



These equations are satisfied only if = 0. The solution is given by the 
following formulae 

X = k'e, + {P{k') + T\k\ k^))e2 + {^{k') + T\k\ A;3))e3, 

(26) 

Here j3 and 7 are arbitrary functions. 

Thus, in the case of the restricted representation (23) for vector a, the 
solution X = x(k), P — const of equations (16) have either of forms (24), 
(25), or (26). 




6. Integration of the equation for vector r 

In all three obtained solutions components and of vector r are 
restricted by the equation 



'2 '2 

J.2 3 
^3 ^3 



1. (27) 



This equation specifies an infinite-dimensional Lie group of area-preserving 
diffeomorphisms (/c^, k^) (r^, r"^). A particular form of the diffeomorphism 
can be obtained by integration of linear differential equation (27) with respect 
to function t^(A;^, k^) for any arbitrary function = T^{k'^, k^) or vice versa. 

Below we use another method of integration of equation (27) in implicit 
form. To this end we perform an incomplete hodograph transformation in 
equation (27). We choose k^ and as new independent variables, and k^, 

as new unknown functions. By the differentiation of identities 

e = K{k\ T\k\ k')) , T^k^ k') = T{k\ T\k\ k')) 

with respect to k^ and k^ we find the transformation of derivatives as 

^2 _ 1 ^2 _ -^fc3 ^3 _ ^3 _ 7fc3X-r2 - Tj.2K^3 

^ K^2 ' ^ K^2 ' ^ ' ^ K^2 

By plugging this formulae into (27) and canceling out the non-zero com- 
mon multiplier Kt-2 we obtain T^s = Kt2. The latter equation is solved by 
introduction of the potential <l>(A;^,r^): 
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Thus, the general solution of equation (27) is implicitly defined by formulae 



(k'y) 



(28) 



with arbitrary smooth function $(A;'^,r^). In particular, the choice 




gives T 



and as 




2A;3sinA;^ = \/2k^ cos k'^ . 



(29) 



Note, that any solution (r^, r^) = (t^(/c^, k^),T^{k'^, k^)) of equation (27) can 
be modified by either of the following ways: 



with arbitrary function G. 

7. Interpretation of obtained solutions 

For the translational surface specified by equation (21) with fixed value 

= c the curve x = cr(A;^,c) will be referred to as the generator, and the 
curve X = t(A;^, c) as the directrix. 

Let us observe solution (24). The generator here is a straight line along 
vector ei, hence contact magnetic surfaces k^ — c are nested cylinders. The 
directrix of cylinders is swept by the vector x = t(A;^, c) whose first com- 
ponent T^{k^,k^) can be chosen arbitrarily, and two remaining components 
should satisfy equation (28). In particular, choice of functions r^, in ac- 
cordance to (29) specifies usual circular cylinders. 

Note, that arbitrary choice of function t^(A;^, k^) does not change contact 
magnetic surfaces. However, the choice of function considerably modifies 
the picture of magnetic lines and streamlines (19) of the flow. In particular, 
by the suitable choice of function plasma flow with fixed contact magnetic 
surfaces can be made either sub-alfvenic (|v| < |B|), alfvenic (|v| = |B|), or 
super-alfvenic (|v| > |B|). 

In the solution (25) the generator is an arbitrary planar curve specified by 
function P{k^). The directrix belongs to a plane spanned by vectors e2 and 



r2 = r2 + G(r3), 



or t'=t' + Gir^) 
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Figure 1: The contact magnetic surface defined by solution (25) with /3 = sinfc^, = 
V2fc3 sinfc^, = V2fc3 cos fc^, fc"^ = 1, F = cos2t^. Continuous and dashed curves on the 
surface are magnetic force fines and streamfines respectively. 



63, and is defined by a solution of equation (27) . Tliis curve can be stretclied 
along ei-direction by a choice of a non-zero function F{t^). Example of the 
contact magnetic surface is given in figure 1. 

Finally, in the solution (26) the generator is a 3D curve specified by 
functions (3 and 7. The directrix is a planar curve given by the solution of 
equation (27). For this solution the modification of a flow is possible due to 
the arbitrariness in the choice of the generator and in the solution of equation 
(27) for the directrix. 

8. Summary 

For incompressible stationary flows of ideal plasma we introduce a con- 
venient curvilinear coordinate system (10), (11), which allows partial inte- 
gration of MHD equations in the form (16) and provides a natural geomet- 
rical description of the flow. For any solution of equations (16) we give 
explicit formulae of streamlines and magnetic force lines (19). We show, that 
Bogoyavlenskij 's infinite-dimensional group of symmetry transformations re- 
duces to a group of scale transformations (17) for the curvilinear coordinates. 
The group preserves contact magnetic surfaces of the flow. Bogoyavlenskij 's 



13 



transformation with discontinuous scale multiplier produces flow with current 
sheet from any continuous flow of plasma. The surface current is explicitly 
given by formula (20). 

By using the curvilinear coordinates we describe stationary flows of ideal 
incompressible plasma with constant total pressure. It is shown, that con- 
tact magnetic surfaces of such flows are translational surfaces, i.e. are swept 
out by translating one curve rigidly along another curve. We give explicit 
examples (24)-(26) of solutions with constant total pressure. The signifi- 
cant functional arbitrariness of solutions allows significant modification of 
described flows. 
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